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BACKSCATTERING  OF  ELECTROMAGNETIC  WAVES  FROM 
A  LAYER  OF  VEGETATION 
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United  States  Army  Engineer  Topographic  Laboratories  ' 

Research  Institute 

Center  for  Theoretical  and  Applied  Physical  Sciences 
Fort  Belvoir,  Virginia  22060 

Abstract 

'  A  theoretical  model  for  the  backscattering  of  electromagnetic  waves 
from  a  layer  of  vegetation  is  computed  using  a  first  order  renormaliza¬ 
tion  technique  to  determine  volume  scattering.  The  vegetation  soil  inter¬ 
face  is  assumed  rough  according  to  the  tangent  plane  approximation  and  the 
scattering  from  this  boundary  is  added  incoherently  to  the  volume  scattering 
result.  The  mean  wave  in  the  vegetation  is  obtained  using  a  bilocal  approxi¬ 
mation  of  the  Dyson's  equation.  A  free  space  dyadic  Green’s  function  is  used 
along  with  a  correlation  function  of  the  dielectric  fluctuations  which  is 

exponential  in  form  and  which  also  possess  different  correlation  lengths 
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&x,  tly,  and  &Z;  in  the  x,  y,  and  z,  directions.  Effective  propagation  constants 
are  obtained  for  both  horizontal  and  vertical  polarizations.  The  scattered 
wave  is  solved  for  using  a  two-dimensional  Fourier  transform  technique  and 
the  boundary  conditions  at  either  end  of  the  vegetation  layer  are  matched. 

t 

The  far  field  backscatter  coefficients  are  computed  for  both  horizontal  and 
vertical  polarizations.  The  mean  and  variance  of  the  dielectric  fluctuations 
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are  calculated  with  the  aid  of  Peake's  model  fcr  the  dielectric  constant  of 
vegetation.  The  theory  is  matched  to  experimental  data  taken  from  a  corn 
field.  The  resulting  values  for  the  correlation  parameters  are  then  used 
to  monitor  the  growth  pattern  of  the  corn  field  over  a  period  of  time. 
Comparison  between  the  theoretical  and  experimental  results  over  this  time 
period  are  shown.  The  theory  is  also  matched  to  experimental  data  from 

spring  and  fall  deciduous  trees.  An  elementary  sensitivity  analysis  is  also 
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shown  in  which  the  change  in  the  }>ackscatter  coefficient  is  plotted  as  a 
function  of  incidence  angle  for  a  given  change  in  one  input  parameter. 
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I.  Introduction 

This  paper  will  present  some  results  of  a  theory  which  was  developed 
for  the  purpose  of  analyzing  the  nature  of  radar  backscattering  from  cer¬ 
tain  types  of  vegetation.  .The  vegetation  is  simulated  by  a  continuous  ran¬ 
dom  medium  and  use  is  made  of  a  first  order  renormalization  technique  to 
calculate  the  radar  backscatter  coefficient.  The  influence  of  an  irregular 
vegetation  soil  interface  is  also  considered  using  a  noncoherent  approach. 
Recently,  Fung1  solved  the  problem  of  scattering  from  a  vegetation  layer 
using  a  scalar  first  order  renormalization  approach.  His  solution,  however, 
did  not  consider  the  existence  of  a  rough  vegetation  soil  boundary.  Fung 
did  consider  an  anisotropic  correlation  function  in  which  the  horizontal 
variation  is  different  from  the  vertical.  Tsang  and  Kong2  solved  the  problem 
of  volume  scattering  from  a  half  space  random  medium  which  contains  lateral 
and  vertical  fluctuations.  A  radiative  transfer  approach  was  used  to  calcu¬ 
late  the  backscattering  cross  sections  up  to  second  order  in  approximation. 
This  enabled  the  cross  polarized  terms  to  be  obtained. 

There  are  two  important  and  somewhat  related  practical  applications 
for  developing  and  analyzing  various  radar  scattering  theories.  The  first 
application  is  radar  image  simulation  of  terrain  features.  In  this  problem 
the  radar  system  parameters  and  terrain  parameters  are  assumed  known  and  used 
to  calculate  a  radar  response  in  the  form  of  a  grey  tone  or  density.  The 
scattering  theories  can  be  used  to  compute  the  radar  backscatter  coefficient 
which  in  turn  is  used  to  calculate  grey  tone.  The  use  of  scattering  theories 
in  this  type  of  application  is  straight-forward  even  though  a  solution  for 
any  one  particular  scattering’ 'problem  may  be  extremely  complicated.  The 
second  application  is  in  the  field  of  remote  sensing  of  terrain  in  which 


i 


the  sensor  response  must  be  utilized  to  determine  various  terrain  parameters. 
The  use  of  the  scattering  theories  for  this  application  is  not  straightfor¬ 
ward.  There  are,  however,  two  important  uses  of  scattering  theories  which 
bear  directly  on  remote  sensing.  The  first  of  these  is  a  parameter  sensitiv¬ 
ity  study.  The  theory  can  be  used  to  analyze  the  influence  of  various  veg¬ 
etation,  terrain  and  radar  parameters  upon  the  sensor  response.  Such  para¬ 
meters  as  surface  roughness,  soil  moisture,  vegetation  height  and  density 
could  be  varied  one  at  a  time  in  order  to  determine  the  influence  on  the 
sensor  response.  This  type  of  analysis  should  lead  to  a  determination  of  what 
radar  parameters  are  most  sensitive  to  certain  terrain  parameter  changes. 

This  type  of  analysis  assumes  the  existence  of  verified  scattering  theories 
which  have  been  developed  and  compared  with  existing  experimental  data. 

The  second  use  of  scattering  theories  would  be  to  analyze  the  radar  response 
for  two  different  types  of  terrain  features  in  an  effort  to  see  if  the  two 
features  could  be  distinguished  from  each  other  on  an  image.  Once  again 
this  would  assume  the  existence  of  verified  scattering  theories  which  have 
been  developed  and  tested  against  experimental  data.  These  applications 
provide  the  incentive  for  developing,  analyzing, and  testing  various  scatter¬ 
ing  theories. 

In  what  follows  we  will  consider  the  geometry  of  the  scattering  problem 
to  be  solved  and  discuss  the  technique  used  for  the  solution.  The  form  of 
the  final  answer  for  the  backscatter  coefficient  will  be  shown,  and  section 
II  will  compare  the  resulting  theory  with  existing  experimental  data.  Also, 
some  results  from  a  sensitivity  of  input  parameters  study  will  be  shown. 

Figure  1  shows  the  scattering  geometry  of  the  vegetation  problem.  A 
plane  wave  with  a  time  harmonic  dependence  of  exp  (jut)  is  incident  from 
free  space  at  an  angle  8^  onto  a  layer  of  vegetation.  The  mean  thickness 


of  Che  vegetation  is  L.  The  vegetation  soil  boundary  is  considered  to  be 
randomly  rough  according  to  the  tangent  plane  approximation.  The  vegetation 
is  simulated  by  a  continuous  random  medium  in  which  e (r)  and  c (r)  represent 
the  three  dimensional  random  dielectric  and  conductivity  fluctuations, 
respectively.  These  fluctuations  consist  of  the  sum  of  an  average  and  a 
fluctuating  component.  The  standard  deviations  of  the  fluctuations  are 
represented  by  nj  and  H2*  The  true  angle  of  refraction  of  the  mean  wave  in 
the  random  medium  is  6e.  The  soil  below  the  vegetation  is  assumed  homogeneous 
with  a  complex  propagation  constant  k.3 .  The  magnetic  permeability  for  all 
three  media  is  assumed  to  be  that  of  free  space.  The  polarization  of  the 
incident  wave  was  taken  to  be  either  horizontal  or  vertical.  A  first  order 
renormalization  method  was  used  to  calculate  the  mean  and  scattered  waves  in 
the  random  medium.  A  solution  was  developed  first  for  the  case  where  the 
vegetation  soil  boundary  is  a  plane  interface.  The  effects  of  an  irregular 
boundary  were  considered  in  a  noncoherent  manner.  The  dielectric  and  con¬ 
ductivity  fluctuation  terms  (e'(r)  and  o'(r))  are  generated  by  statistically 
homogeneous  random  processes.  The  mean  of  each  process  is  zero  and  the 
correlation  functions  are  identical  and  exponential  in  form  with  1 Zy  and  Zz 
representing  the  correlation  distances  in  x,  y,  and  z.  The  correlation 
functions  have  been  chosen  to  be  anisotropic.  It  is  believed  that  this 
representation  with  unequal  correlation  distances,  is  closer  to  reality  than 
an  isotropic  correlation  function.  This  is  because  the  size  of  the  vegetation 
scatterers  in  a  horizontal  plane  is  generally  not  the  same  as  the  size  of  the 
scatterers  in  a  vertical  plane.  The  mean  wave  in  the  random  medium  is  determined 
'from  the  bilocal  approximation' of  the  Dyson's  equation.  Plane  wave  solutions 
to  the  Dyson’s  equation  were  obtained  using  the  form  of  a  free  space  dyadic 
Green's  function.  The  free  space  propagation  constant  was  replaced  with  the 
average  background  propagation  constant  of  the  random  medium.  A  three 


dimensional  Fourier  transform  of  the  Green's  function  was  used.  A  solution 
for  the  2  component  of  the  effective  propagation  constant  of  the  mean  wave 
was  obtained  for  both  horizontal  and  vertical  polarizations.  Once  the  mean 
wave  had  been  calculated  and  the  appropriate  boundary  conditions  matched  for 
the  mean  waves  in  all  three  media,  then  the  scattered  wave  in  the  random 
medium  was  computed  using  the  first  order  renormalized  equation  for  the 
scattered  wave.  This  scattered  wave  was  computed  using  a  two  dimensional 
Fourier  transform  technique.  This  also  allowed  the  calculation  of  the 
scattered  waves  in  air.  The  necessary  boundary  conditions  were  matched  and 
the  backscatter  coefficient  was  computed  for  horizontal  and  vertical  polar¬ 
izations.  The  influence  of  the  rough  boundary  between  the  vegetation  and  the 
soil  was  considered  apart  from  the  volume  scattering  solution,  using  the 
tangent  plane  method.  The  backscatter  coefficient  for  the  rough  surface 
scattering  had  to  be  modified  by  the  attenuation  through  the  vegetation. 

This  result  was  then  added  to  the  volume  scattering  solution  to  obtain  a 
final  answer  for  the  backscatter  coefficient.  The  form  of  this  result  is 
given  below: 

°HH  *  °HHS  exP(-4oelLsec0el)'K7HHV 

°VV  =  °VVS  exP(~^ae2Lsec6e2)+0vVv  (2) 

where: 

o§hs  an<^  a$VS  are  t*ie  rough  surface  scattering  coefficients  for 
horizontal  and  vertical  polarizations. 

t 

ae^  and  ae2  are  the  attenuation  portions  of  the  effective  propagation 

constant  for  horizontal  and  vertical  polarizations,  respectively. 

6e^  and  0e2  are  the  true  angles  of  refraction  of  the  mean  wave  in 
the  random  medium  for  horizontal  and  vertical  polarizations, 
respectively. 
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°HHv  an<*  0Wv  are  t*ie  vo^ume  scattering  coefficients  for  horizontal 
and  vertical  polarizations. 

An  elementary  permittivity  model  was  developed  to  relate  certain  para¬ 
meters  of  the  random  medium  to  the  parameters  of  actual  vegetation. 

Use  was  made  of  a  model  provided  by  Peake  and  Oliver3  for. the  complex 
dielectric  constant  of  vegetation.  This  in  turn  was  used  to  help  compute 
values  for  the  means  and  variances  of  the  dielectric  and  conductivity 
fluctuations. 

II.  Numerical  Results 

This  section  will  show  the  results  of  sample  calculations  for  equations 
(1)  and  (2)  and  will  also  show  the  results  of  some  sensitivity  calculations 
where  the  change  in  backscatter  coefficient  is  analyzed  for  a  given  change 
in  one  input  parameter.  The  input  parameters  for  the  model  are: 

1.  Angle  of  incidence  (ei) 

2.  Fraction  of  water  by  weight  in  the  vegetation  (F) 

3.  Volume  of  vegetation  divided  by  the  total  volume  of 

vegetation  plus  air  (Ry) 

4.  Correlation  distance  in  the  x  direction 

5.  Correlation  distance  in  the  y  direction  (£y) 

6.  Correlation  distance  in  the  z  direction  (£z) 

7.  Mean  thickness  of  the  vegetation  layer  (L) 

8.  Relative  dielectric  constant  of  the  soil 

below  the  vegetation  (Eg) 

9.  Conductivity  of  the  soil  below  the  vegetation  (o^) 

10.  Frequency  , '  (f ) 

11.  Ratio  of  the  standard  deviation  of  the  rough 

surface  fluctuations  to  the  correlation  distance 
of  the  fluctuations  (ms) 
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The  output  of  the  computer  calculations  is  the  backscatter  coefficient  in 
decibels.  Figures  2  through  7  match  the  theory  with  experimental  data  taken 
from  corn  at  three  different  times  during  the  growing  season.  The  peak 
which  occurs  at  normal  Incidence  is  due  to  rough  surface  scattering  and 
changes  according  to  the  thickness  of  the  layer  and  the  vegetation  moisture 
content.  Figures  8  and  9  match  the  theory  to  experimental  data  taken  from 
spring  and  fall  deciduous  trees.  These  results  are  also  compared  to  some 
results  derived  for  corn.  Figure  10  presents  a  study  of  a°  variations  with 
the  fraction  of  water  by  weight  in  the  vegetation.  For  small  angles  of 
incidence  (less  than  10°)  increasing  F  results  in  a' slightly  lower  value  of 
o°.  This  is  due  to  the  fact  that  surface  scattering  is  dominant  at  these 
angles  and  Increases  in  F  yield  more  attenuation  for  the  mean  wave.  At 
larger  angles  of  incidence  (greater  than  20°)  volume  scattering  is  dominant 
and  increases  in  F  increase  the  level  of  o°.  Figure  11  presents  a  study  of 
o°  variations  with  L.  Figure  12  shows  a  study  of  o°  variations  with  £x*  For 
angles  of  incidence  less  than  approximate  50°,  an  increase  in  ■£*  results  in 
larger  values  of  o° .  For  angles  of  incidence  greater  than  50°  an  increase 
in  -£x  results  in  a°  falling  off  faster  with  incidence  angle.  For  angles  of 
incidence  less  that  10°  changes  in  Zx  have  no  influence  on  a°  since  the 
surface  scattering  is  dominant.  Figure  13  presents  a  study  of  a°  variations 
with  Zz.  Increasing  Zz  results  in  increasing  the  level  of  o°  except  for 
angles  of  incidence  less  that  10°.  Figure  14  shows  a  study  of  o°  variations 
with  Eg.  In  this  figure  the  change  in  o°(Aa°)  is  plotted  as  a  function  of 
incidence  angle  for  a  change  in  Eg  of  40.  It  can  be  seen  that  for  small 
angles  of  incidence  this  results  in  large  changes  in  o°.  As  the  incidence 


angle  Increases,  however,  Ao°  rapidly  decays  toward  zero  as  the  lower  surface 
has  less  influence  on  the  result.  Figure  15  shows  a  study  of  o°  variations 
with  £x*  The  change  in  o°  is  plotted  as  a  function  of  incidence  angle  for 
a  change  in  t equal  to  one  centimeter.  For  angles  of  incidence  less  than 
20°,  Ao°  is  practically  zero.  At  these  angles  the  rough  surface  scattering 
is  most  inportant  and  changes  in  -£x  have  little  or  no  affect  on  o°.  As  the 
angle  of  incidence  increases;  Aa°  swings  downward  to  a  value  of  approximately 
-4db  at  80°. 

We  will  conclude  this  section  with  a  brief  discussion  of  the  limitations 
and  difficulties  with  the  developed  theory.  It  appears  valid  to  simulate 
a  region  of  vegetation  with  a  continuous  random  medium  although  it  is  not 
certain  as  to  how  well  the  first  order  renormalization  technique  does  in 
solving  the  problem.  It  is  not  cleat  how  much  multiple  scattering  is  being 
considered  and  it  is  not  even  clear  as  to  how  much  multiple  scattering  must 
be  considered.  A  free  space  dyadic  Green's  function  was  used  in  solving 
the  Dyson's  equation  and  what  should  have  been  used  was  a  Green's  function 
applicable  to  a  layered  problem.  An  anisotropic  correlation  function  was  used, 
however,  this  did  not  result  in  a  depolarization  term.  The  existence  of 
depolarization  is  clearly  evident  from  the  experimental  data.  The  reason 
for  this  depolarization  is  as  yet  unknown.  A  depolarization  term  could  be 
obtained  by  computing  the  scattered  field  to  a  second  order  approximation.  It 
could  also  be  obtained  by  initially  allowing  for  an  anisotropic  random  medium. 
Which  approach  is  correct  is  unclear  at  this  time. 
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III.  Conclusions 


The  following  conclusions  have  been  made  as  a  result  of  this  work: 

1.  A  theory  has  been  developed  to  explain  electromagnetic  wave 
scattering  from  a  vegetation  layer  possessing  an  irregular  vegetation 
soil  boundary. 

2.  The  theory  developed  was  for  like  polarized  (HH  or  W)  components 
only  and  no  result  was  obtained  for  the  depolarized  components.. 

3.  For  certain  types  of  vegetation  such  as  corn  it  was  found  that 
the  irregular  vegetation  soil  boundary  dominates  the  backscattering 
result  for  angles  of  incidence  between  0°  and  20°. 

4.  In  general  the  affect  of  the  rough  surface  boundary  between  the 
vegetation  and  soil  increases  with  decreases  in  frequency,  vegetation 
moisture  content,  vegetation  volume,  and  layer  thickness. 

5.  The  o°  versus  incidence  angle  curve  was  shown  to  be  sensitive  to 
very  slight  changes  in  the  correlation  distance  in  z. 

6.  The  predictability  of  the  0°  versus  incidence  angle  curve  depends 
upon  a  very  detailed  knowledge  of  the  statistical  properties  of  the 
dielectric  fluctuations  of  the  vegetation  and  surface  roughness 
properties  of  the  soil  below.  Such  knowledge  as  this  for  particular 
vegetation  features  does  not  exist  at  the  present  time.  Future  work 
should  include  attempts  to  obtain  this  detailed  understanding  if 
theoretical  models  are  to  have  ultimate  usefulness  in  predicting 
scattering  from  vegetation  features. 
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RNGLE  DF  INCIDENCE  IN  DEGREES 


Polarization:  Horizontal 
F  *  0.815  f  =  8.6GHz 
R„  =  0.0002 
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zation:  Vertical 
815  f  =  8.6GH 
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RNELE  DF  INCIDENCE  IN  DEGREES 


Polarization:  Horizontal 
F  =  0.529  f  =  8.6GHz 
Rv  =  0.0002 

lx  =  ly  -  5  m  Iz  -  0.095 


ANGLE  OF  INCIDENCE  IN  DEGREES 


Vertical 
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Polarization:  Horizontal  f  =  8.6GHz  Corn  Ground  Truth 

F  =  0.895  R„ 
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Polarization:  Horizontal 


ANGLE  DF  INCIDENCE  IN  DEGREES 


Polarization:  Horizontal 
F  =  0.8  f  =  8GHz 

Rw  =  0.0002 


RNGLE  DF  INCIDENCE  IN  DEGREES 


Polarization:  Horizontal 
F  «  0.82  f  «  8.6GHz 
Rw  =  0.0008 


